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Energies of (2s 2 2p 5 3Z3/ ; ) and (2s2p 6 3l3l') states for sodiumlike ions with Z =14-100 are evaluated 
to second order in relativistic many-body perturbation theory starting from a neonlike Dirac-Fock 
potential. Second-order Coulomb and Breit interactions are included. Correction for the frequency- 
dependence of the Breit interaction is taken into account in lowest order. The Lamb shift correction 
to energies is also included in lowest order. Intrinsic particle-particle-hole contributions to energies 
are found to be 20-30% of the sum of one- and two-body contributions. Comparisons are made 
with available experimental data. These calculations are presented as a theoretical benchmark for 
comparison with experiment and theory. 

PACS numbers: 31.15.Ar, 31.15.Md, 31.25.Jf, 32.30.Rj 
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I. INTRODUCTION 



Core-excited states of Na-like ions have two valence 
electrons outside and one hole inside a closed n=2 core 
and, therefore, present an excellent model for studying 
strong correlations in systems with closely-spaced levels. 
We carry out second-order relativistic many-body per- 
turbation theory (RMBPT) calculations for Na-like ions 
starting from a ls 2 2s 2 2p 6 Dirac-Fock (DF) potential. All 
possible 21 holes and 3l'3l" particles leading to 121 odd- 
parity and 116 even-parity 3l'3l"[Ji]2l~ 1 ( J) states are 
considered. We calculate energies of the 237 core-excited 
states together with the 31 singly-excited states in Na-like 
ions with nuclear charges Z ranging from 14 to 100. 

Many experimental values for the energy levels and 
fine-structure intervals are now available up to very high 
nuclear charge (Z=57) for 3l'3l"[J 1 ]2l- 1 {J) levels; ad- 
ditionally, experimental rates for some transitions be- 
tween these levels are available. The objective of this 
paper is to present a comprehensive set of calculations 
for 3l'3l"[Ji]2l~ 1 (J) energies to compare with previous 
calculations and experiments for the entire Na isoelec- 
tronic sequence. Most earlier measurements and calcu- 
lations focused on the 3s 2 [ 1 5]2p _1 2 Pj states and low- 
lying 3s3p[ 3 P}2p- 1 4 Dj and 3s3d[ 3 D}2p- 1 4 Fj levels. 
Very few results exist for other 3p 2 2p~ 1 , 3S3P2P" 1 , and 
3s3d2p~ 1 levels. The large number of possible transitions 
have made experimental identification difficult. Experi- 
mental verification should become simpler and more re- 
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liable using this more accurate and complete set of cal- 
culations. 

We start with a brief overview of previous theoreti- 
cal and experimental studies of properties of core-excited 
states in sodium-like ions. Energies, collision strengths, 
and oscillator strengths for 2s 2 2p 6 3l - 2s 2 2p 5 3l'3l" and 
2s 2 2p 6 3l - 2s2p 6 3l'3l" transitions in Na-like ions with 
nuclear charges 22 < Z < 62 were calculated by Zhang 
et al. 0, where energies and mixing coefficients were 
obtained using the Cowan code Q with a scaling fac- 
tor (Z - 3.9) /Z. In the paper by Chen §, Auger 
and radiative transition energies and rates were calcu- 
lated for 18 ions with nuclear charges in the range 18 < 
Z < 92 using the multiconfiguration Dirac-Fock (MCDF) 
method. Wavelengths, radiative transition rates, Auger 
rates, satellite intensity factors, and fluorescence yields 
were presented by Nilsen B for dielectronic satellites 
of 14 neon-like ions (Ar 8 + - W 64+ ) including 2s 2 2p 6 3l- 
2s 2 2p 5 3l'3l" and 2s 2 2p 6 3l -2s2p 6 3l'3l" transitions. The 
calculations in Q were made using the YODA code which 
is based on multiconfiguration relativistic functions for 
the bound states and distorted-wave Dirac wave func- 
tions for continuum electrons. New comprehensive theo- 
retical analyses of the 2s 2 2p 5 3lnl' and 2s2p 6 3lnl' (n=3-8) 
doubly-excited states in Fe 15+ and Cu 18+ were presented 
by Bruch et al. ||, g|. In these papers, energy levels, ra- 
diative transitions probabilities, and autoionization rates 
were calculated by using the Cowan and YODA codes. Re- 
cently, R-matrix calculations of electron-impact collision 
strengths for excitation from the inner L-shell into dou- 
bly excited states of Fe 15+ were presented by Bautista . 
Energy positions of the 2s 2 2p 6 3l and 2s 2 2p 5 3l3l' states 
were calculated in that work using the superstructure 
code of Eissner et al. ||. It was shown in 0] that dis- 
agreement between calculations and data recommended 
by Sugar and Corliss |j| and Shirai et al. |l(J ranges from 
0.5% to 5%. 
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Studies of the 2s 2 2p 5 3lnl' configurations in Na-like ions 
are of continuing interest from both theoretical and ex- 
perimental points of view. Experimentally, these con- 
figurations are studied by photon and electron emission 
spectroscopy. To our knowledge, the first measurements 
of 2p~ 3s and 2p—3d transitions in Na-like Fe (which in- 
cluded the 2s 2 2p 5 3l3l' configurations) using photon emis- 
sion spectroscopy were done by Burkhalter et al. [Ill in 



1979. The soft x-ray spectra produced by a focused laser 



source and a vacuum spark device were obtained in |ll] 
using a grazing incidence spectrograph. Identification of 
the spectra was based on theoretical calculations using 
the Cowan code. The 2s 2 2p 5 3s3p 4 D 7/2 - 2s 2 2p 5 3s3d 
4 Gg/2 transition in FeXVI was observed by the beam-foil 
technique by Jupen et al. jl2|]. The identification of this 
transition was based on a study of the difference in wave 
numbers between experimental values and theoretical 
predictions obtained from the Cowan code in Na-like ions 
from S VI to Cu XIX. Auger electron spectra of Na-like Fe 
ions excited in collisions of 170-keV Fe XVIII on He and 
Ne were studied by Schneider et al. |l3| . The dominant 
spectral structures were due to Auger decay of states with 
2s 2 2p 5 3l3l' configurations in FeXVI. Special attention 
was paid to the metastable level 2s 2 2p 5 3s3p A D 7 / 2 which 
was studied in detail; the absolute energy of this level was 
measured. Three of the 17 lines observed in this experi- 
ment were identified with states from 2s 2 2p 5 3lAl' config- 
urations. This analysis was based on theoretical predic- 
tions using the multiconfigurational Dirac-Fock method 
including the transverse Breit interaction. A few years 
later this analysis was extended to include an isoelec- 
tronic study of double electron capture in slow ion-atom 
collisions of F-like ions (Si 5+ , Ar 9+ , Sc 12+ , Ti 13 +, Fe 17+ , 
and Cu 20+ ) with He atoms [0. It was reported that 
the metastable level 2s 2 2p 5 3s3p A D 7 / 2 was observed in 
all Auger spectra except Si 3+ . Recently in |15|j , a study 
of the double electron capture in low-energy Fe 17+ + He 
collisions was conducted. Except for the confirmation of 
the identification of the 12 peaks in the 225 - 340 eV 
interval, predictions were made for identification of four 
additional peaks in the 585 - 670 eV interval. A com- 
parison of solar and tokamak FeXVII x-ray spectra with 
synthetic spectra including dielectronic satellites with 
2s 2 2p 6 nl - 2s 2 2p 5 3snl and 2s 2 2p 6 nl - 2s 2 2p 5 3dnl tran- 
sitions was presented by Phillips et al. Q. Atomic data 
were calculated in Phillips et al. [^6| using the super- 
structure code. Recently, the three Fe 15+ inner-shell 
satellite lines were thoroughly studied at the Lawrence 
Livermore National Laboratory electron beam ion trap 
ebit-ii by Brown et al. Q . 

In the present paper, we use RMBPT to determine en- 
ergies of 2s 2 2p 5 3l3l' and 2s2p 6 3l3l' states for Na-like ions 
with nuclear charges in the range of Z = 14-100. Our 
calculations are carried out to second order in perturba- 
tion theory and include second-order Coulomb and Breit 
interactions. Corrections for the frequency-dependent 
Breit interaction are taken into account in the lowest or- 
der. Screened self-energy and vacuum polarization data 



given by Blundell |18 are used to determine the QED 
correction. 

Our perturbation theory calculations are carried out 
using single-particle orbitals calculated in the DF poten- 
tial of the Ne-like core. As a first step, we determine and 
store the single-particle contributions to the energies for 
three n—2 hole states (2s, 2p 1 / 2 , and 2p 3 / 2 ) and the five 
n=3 valence states (3s, 3pi/ 2 , 3p 3 / 2 , 3d 3 / 2 ^ and 3^5/2) in 
lowest, first, and second orders. Next, we evaluate and 
store the 155 two-particle (3131' J\H eS \3l"3l"' J) matrix 
elements and the 178 hole-particle (2131' J\H cS \2l"3l"' J) 
matrix elements of the effective Hamiltonian in the first 
and second orders. It should be noted that these one- 
particle, two-particle, and hole-particle matrix elements 
were used previously to evaluate energies of the 3131' lev- 
els in Mg-like ions [|19| and energies of the 2l~ x 3l' lev- 
els in neonlike ions pOR , Finally, second-order particle- 
particle-hole matrix elements are evaluated. Combining 
these data using the method described below, we cal- 
culate one-, two-, and three-body contributions to the 
energies of Na-like ions. 

The present calculations are compared with experi- 
mental and predicted results from Refs. [O, \\Jl BU E2). 



II. METHOD 

The RMBPT formalism developed previously ||| ||, 
p5| for B-like and Al-like ions is used here to describe the 
perturbed wave functions and to obtain the second-order 
energies [^3| of the core-excited states of Na-like ions. 
The particle-particle-hole contribution to the energies of 
those states is, however, different from the three-particle 
contribution to the energies of B-like and Al-like ions, 
and is derived in the Appendix. Including hole-particle 
matrix elements in Na-like ions leads to 237 313V21" 1 
states instead of 148 313V31" states of Al-like ions or 15 
212V21" of B-like ions and, consequently, to more labori- 
ous numerical calculations. 



A. Model space 

The model space for core-excited 3l3l'2l~ 1 states of 
Na-like ions includes 121 odd-parity states consisting of 
25 J=l/2 states, 37 J=3/2 states, 31 J=5/2 states, 19 
J=7/2 states, 7 J=9/2 states, and two J=ll/2 states. 
Additionally, there are 116 even-parity states consisting 
of 25 J=l/2 states, 35 J=3/2 states, 31 J=5/2 states, 
17 J=7/2 states, 7 J=9/2 states, and one J=ll/2 state. 

The evaluation of the second-order energies for the 
3l3l'[Ji]2l~ 1 ( J) states in Na-like ions follows the pat- 
tern of the corresponding calculation for Mg-like and Ne- 
like ions given in Refs. Jl^, ^0). In particular, we use 
the second-order one-particle and two-particle matrix ele- 
ments for Mg-like ions calculated in 1 19 and hole-particle 
matrix elements for Ne-like ions calculated in [Eel . These 
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matrix elements are recoupled as described below to ob- 
tain the one-body and two-body contributions for Na- 
like ions. We will discuss how these matrix elements are 
combined to obtain the one-body and two-body contri- 
butions to energies of Na-like ions. We refer the reader 
to Refs. (TH " 



20| for a discussion of the evaluation of 
the one-body and two-body matrix elements. An intrin- 
sic particle-particle-hole diagram also contributes to the 
second-order energy for Na-like ions. The expression for 
particle-particle-hole diagram differs from the expression 
for the three-particle diagram given in Ref. |23| owing to 
different angular parts. The expression for the intrinsic 



order contribution S' ' can be about half of B {2) . Fi- 
nally, we observe some small deviations from smooth in- 
dependencies of E^ a nd B^ in Figs. ||a and ||b. As we 



particle-particle-hole diagram is derived in Appendix. 



can see from Eq. ( A10[ ) in the Appendix, the denominator 
e(v) + e(w) — e(a) — e{n) in the expression for the three- 
body diagram contains two valence energies s(v), s(w), 
one hole energy e(a), and one virtual state energy e[n). 
The sum e(3s) + e(3s) — e(2p) is positive; therefore, the 
denominator can be close to zero when e{n) is positive. 
In such a case, a singularity arises from the contribution 
of the continuous part of spectra to the sum over states 

(2) (2) 

for the E$ and matrix elements, resulting in sharp 
features in the corresponding Z-dependencies. 



B. Energy-matrix elements 
In Fig. pla, we present the one-body, E^\a,a), two- 

(2) (2) 

body, i?2 (a, a), and three-body, E$ (a,a) Coulomb 
contributions to the diagonal matrix element with a = 
3s 1 / 2 3si/2[0]2p 1 / 2 (1/2) as functions of Z. Here and in 
the sequel, we drop the superscript —1 on the n = 2 
hole state for simplicity. We see from Fig. pla that the 

(21 

value of £3 (a, a) is almost constant for all Z. The 

one-body contribution is the largest one. The sum 

(2) (2) 
of E 2 (a, a) and E\ (a, a) values is almost equal to 

zero for small-Z ions; however, the three-body contri- 
bution becomes more important with increasing Z. As 
a result, the total second-order Coulomb contribution 

E^(a,a)=E[ 2 \a,a)+E^ ) (a,a)+E^\a,a) differs from 
the one-body contribution and is almost constant for 

all Z. In Fig. [lib, we present the one-body, 5^ (a, a), 

(2) (2) 
two-body, B 2 '(a, a), and three-body, B\ (a, a) Breit- 

Coulomb contributions to the diagonal matrix element 
with a — 3s 1 / 2 3s 1 /2[0]2p 1 /2 (1/2) as function of Z. 
The one-body contribution is the largest among other 
contributions and the total Breit second-order contri- 
bution B^{a,a)^B {2 \a,a)+B^ 2 \a,a)+B {2 \a,a) is al- 
most equal to the one-body contribution. The one-body 
contribution vanishes for non-diagonal matrix elements. 

In Fig. ga, we illustrate the contributions E^ 

(2^ C2) 

and E% for two non-diagonal elements E\ (a, 6) 

( 2) 

and El (b,a) with a=3s 1 / 2 3s 1 / 2 [0]2p 1 / 2 (l/2) and 
&=3pi/ 2 3pi/ 2 [0]2p!/ 2 (l/2). The corresponding Breit con- 
tributions are plotted in Fig. ||b. We see from Fig. ||a 
that the ratio of E {2) to E {2) is about 30-50% for high-Z 
ions. Fig. ||a illustrates the asymmetry of the nondiag- 
onal energy matrix elements for low-Z ions. The asym- 
metry of the two-particle second-order matrix elements 
in RMBPT calculation was discussed in Ref. [2(J. Val- 

(2) (2) 

ues of E 2 ' and E^ ' have different signs and the total 

second-order contribution E^ is about half of E 2 \ The 
behavior of non-diagonal second-order Breit contribution 
is shown in Fig. ||b. The B 2 and B^ increase with Z; 
however, they have opposite signs. The total second- 



er. Eigenvalues and eigenvectors for core-excited 
states 

After evaluating the energy matrices, we calculate 
eigenvalues and eigenvectors for states with given val- 
ues of J and parity. There are two possible methods to 
carry out the diagonalization: (a) diagonalize the sum 
of zeroth- and first-order matrices, then calculate the 
second-order contributions using the resulting eigenvec- 
tors or (b) diagonalize the sum of the zeroth-, first-, and 
second-order matrices together. Following Ref. p3| ), we 
choose the second method here. 

The importance of second-order contributions to the 
energies is illustrated in Fig. g. In this figure, the second- 
order energy E^ 2 \ the first- and second-order Breit con- 
tributions B^ and B^ 2 \ and the QED contribution 
-E"lamb are plotted as functions of nuclear charge Z for 
the 3s 2 { 1 S)2p 2 P 1/2 (a) and 3s3p{ 3 P)2p A D 7/2 (b) states 
of Na-like ions. We can see from the figure that E^ 
dominates up to Z=34 and Z—37 for the 3s 2 2p 2 Pi/ 2 
and 3s3p( 3 P)2p A D 7 / 2 states, respectively. The QED 
contribution -Elamb is smaller than the E^ energy up 
to Z=50 for the energy level shown in Fig. ga and Z=60 
for the energy level shown in Fig. gb. The second-order 
Breit energy B^ is smaller than the second-order energy 
E^ up to Z=97 (Fig. ga) and Z=83 (Fig. gb), respec- 
tively. 

The sharp features in the curves shown in the 
Fig. ga are explained by strong mixing of states in- 
side of the odd-parity complex with J=l/2. The main 
contribution to the 3s 2 ( 1 S)2p 2 Pi/ 2 level comes from 
the 3si/ 2 3si/ 2 [0]2p!/ 2 state for ions with small Z up 
to Z=A8. Starting from Z=49, the main contribu- 
tion to this level comes from 3p!/ 2 3p3/ 2 [l]2p 3 / 2 and 
3pi/ 2 3p 3 / 2 [2]2p 3 / 2 states instead of 3si/ 2 3s 1 / 2 [0]2p 1 / 2 
state. 

The particle-particle-hole levels studied in this work 
may be divided to two groups by their energies: those 
with energies in the range of 15-23 a.u and those 
with energies greater than 25 a.u. The first group 
of levels includes 3s 1 / 2 3s 1 / 2 [0]2p J -, 3s 1 / 2 3c?j[Ji]2p J v, 
3pj3pj>[Ji]2pj" , and 3si/ 2 3p J [Ji]2p J / levels (64 levels to- 
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FIG. 1: Second-order Coulomb (a) and Breit (b) energy contributions to diagonal matrix elements in Na-like ions. 
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FIG. 2: Second-order Coulomb (a) and Breit (b) energy contributions to non-diagonal matrix elements in Na-like ions. 



tal). The second group includes the remaining 173 lev- 
els. The first group of levels has been studied experimen- 
tally, while there are no experimental data for the second 
group. 

When starting calculations from relativistic DF wave 
functions, it is natural to use jj designations for uncou- 
pled energy matrix elements; however, neither jj nor 
LS coupling describes physical states properly, except 
for the single-configuration state 3p3/23d,5/2[ / t]2p 3 /2 = 
3p3d[ 3 F]2p 3 Gn/2- Both designations are given in our 
subsequent tables. 

The Z-dependence of the eigenvalues of the 
3lj3l'j'[J 1 ]2l"j"(J) odd-parity states with J=l/2 
is given in Figs. ||a, b, c, and d, where the energies are 
divided by (Z — 3) 2 . Scaled energies E/(Z — 3) 2 decrease 



with Z for the first group of levels (Figs. [|a, b) and have 
parabolic Z dependence for the second group (Figs. ||c, 
d). It is evident from Figs. Ha, b that there is a strong 
mixing between 3s 1 / 2 3s 1 /2[0]2p 1 /2, 3s 1 / 2 3c?j[Ji]2p ; ,/, 
and 3pj3pj'[Ji]2pj>> levels. We already mentioned that 
strong mixing happens between 3si/23si/2[0]2p 3 / 2 , 
3pi/23p 3 /2[l]2p 3 /2 and 3pi/ 2 3p3/2[2]2p 3 /2 states for 
Z=48-49. 

It should be noted that the LS designations are chosen 
by comparing results with available experimental data 
for low-Z ions. We used the labeling of levels from 
Refs. p, p), where the identification of levels in Fc 15+ 
and Cir 8 " was based on the Cowan code. We already 
mentioned that there are no measurements for the second 
group of levels, therefore, we do not give LS designations 
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FIG. 3: Contributions to the energies of the 3s 2 ( 1 S)2p 2 Pi/2 (a) and 3s3p( 3 P)2p 4 D 7 / 2 (b) levels of Na-like ions. 
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FIG. 4: Energies (E/(Z — 3) 2 in cm x ) of odd-parity states with J— 1/2 as functions of Z. 
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for the levels shown in Figs. |]c:, d. The change of rela- 
tive magnitude of mixing coefficients is illustrated by the 
behavior of the six levels in Fig. ||a and the five levels 
in Fig. [IJd. The crossing of levels inside of a complex 
(set of levels with the same J and parity) is forbidden by 
Wigner-Neumann theorem p~7f . 



III. COMPARISON OF RESULTS WITH 
OTHER THEORY AND EXPERIMENT 



Identification of spectral lines originating from the tran- 
sitions between core-excited configurations in Na isoelec- 
tronic sequence was made from analyses of beam-foil 
spectra in Jl^| . Wavelengths A cxp t from are shown in 
Table [Fv| together with our RMBPT data Armbpt- We 
also present the predicted data Agt from Ref. (l2) in the 
same table. The predicted wavelengths were obtained in 
|]l2| from the fit of the difference between experimental 
data and Cowan code calculations. Our ah initio calcu- 
lations are in good agreement with those predictions. 



We have calculated energies of the 121 odd-parity 
states and the 116 even-parity excited states for Na- 
like ions with nuclear charges ranging from Z=1A to 
100. In Table |, we illustrate our theoretical results giv- 
ing energies of the 38 odd-parity states 3s 2 ( 1 S)2p 2 Pj, 
3s3d( 13 D)2p 2 - 4 L,j, and 3p 2 { 1 ' 3 L')2p 2 ' 4 L. T and the 18 
even-parity states 3s3p( 13 P)2p 2,4 Lj for eight Na-like 
ions with nuclear charges ranging from Z=26 to Z=57. 
We also include the ground state 3s 2 Si/ 2 energies in the 
table. 



A. Excitation energies 

Comparisons of our RMBPT energies with other the- 
oretical and experimental data are too voluminous to 
include here; therefore, we present several examples of 
comparisons for selected levels and ions. In Table III our 
theoretical results for 3s 2 { 1 S)2p 2 P J: 3s3p{ 1 ' 3 P)2p^ 4 L J , 
3p 2 { 1 ' 3 L')2p 2A Lj, and 3s3d{ 1 ' 3 D)2p 2 > 4 Lj levels of Na- 
like Fe are compared with predicted Auger energies given 
by Bliman et al. 15 1. It should be noted that the pre- 
dicted Auger energies in Ref. |H| differ by 1 eV from ex- 
perimental measurements presented more than 15 years 
ago by Schneider et al. 13 , where the method of zero- 
degree Auger spectroscopy was used to study the excita- 
tion and decay of Auger states formed in 170-keV Fe 17+ 
ions with He and Ne target atoms. Our RMBPT energies 
are in much better agreement with recent results from 
Jl5[ than with experimental measurements from ]13|] . 

In Table |in[ our RMBPT wavelengths for Na-like Ag 
and La are compared with experimental results presented 
by Beiersdorfer et al. |2l], p2| . Spectra of the n=3 to n—2 
transitions in neonlike silver obtained from the Princeton 
Large Torus were recorded with high-resolution Bragg- 
crystal spectrometer. The measurements covered the 
wavelength interval 3.3-4.1 A. Eight lines of Na-like Ag 
were observed. As one can see from Table III, our theo- 



B. 



Fine structure of the 2 L and 4 L terms in 
core-excited states of Na-like ions 



retical calculations are in excellent agreement with re- 
sults from Ref. [^J. Measurement of a single line in 
Na-like La was carried out with high-resolution crystal 
spectrometer on the PLT tokamak by Beiersdorfer et al. 
p2| . Our RMBPT result agrees with experimental wave- 
lengths from pa to five significant figures. We list our 



theoretical data for other lines in La 46+ in Table III. 

The 3s3p{ 3 P)2p 4 D 7/2 -3s3d( 3 D)2p 4 F 9/2 transitions 
in P 4+ -Zn 19+ ions were investigated by Jupen et al. [p"2"[. 



No measurements of fine-structure intervals were made 
by observing the wavelength differences between transi- 
tions within the doublet or quartet states. The intervals 
are obtained by subtracting the excitation energies of the 
corresponding levels (see, for example, Auger data in Ta- 
ble [n]) . These fine structures are quite regular through- 
out the isoelectronic sequence, as seen from Fig. |^a. In 
this figure, we present the fine-structure splitting di- 
vided by (Z - 3) 2 for the doublet term 3s 2 { 1 S)2p 2 P. 
Small deviation from smooth line of Fig. |^a in the re- 
gion of Z=48-49 is due to the strong mixing in the odd- 
parity complex with J=l/2 as illustrated in Fig. ||a. 
The 3s 2 ( 1 S)2p 2 P^/ 2 level has smaller energy than the 
3s 2 ( 1 S)2p 2 P\/2 level in contradiction with Lande in- 
terval rules. The same situation occurs for other fine- 
structure intervals in Figs. [|b, c, and d. The 4 P 3/2 level 
has the smallest energy among 3p 2 ( 3 P)2p 4 Pj levels, the 
next level is 4 P5/2, and the 4 -Pi/2 level has the largest 
energy. The value of {Ei/ 2 — E 3 / 2 )/(Z — 3) 2 decreases 
with Z while the value of (E 5 / 2 — E 3 / 2 ) / (Z — 3) 2 increases 
with Z. As a result, energies of 4 -Pi/2 and 4 P§/ 2 become 
almost equal for high Z . The energy difference between 
the 4 -D 3 / 2 and 4 D 5 / 2 levels in the 3s3p( 3 P)2p 4 Dj quar- 
tet term is very small, as shown in Fig. ||c. We find that 
the 4 Dj quartet terms are ordered as follows: 4 D 1 / 2 , 
A D 7 / 2 , 4 £>3/2, and 4 -D5/2- Energy differences between 
the 4 F 5/2 , 4 F 7/2 , and 4 F 9/2 levels in the 3s3d{ 3 D)2p 4 Fj 
quartet term are very small, as shown in Fig. pp. The 
unusual splitting is caused by mixing of states within the 
complexes with different J. Further experimental confir- 
mation would be very helpful in verifying the correctness 
of these occasionally sensitive mixing parameters. 



IV. CONCLUSION 

A systematic second-order RMBPT study of the en- 
ergies of the core-excited states of Na-like ions has been 
presented. RMBPT gives results in good agreement with 
experimental and predicted data. It would be beneficial if 
experimental data for other highly-charged Na-like ions 
were available. At the present time, there are almost 
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TABLE I: Energies of Na-like ions relative to the ground state of Ne-like ions in cm for ions with Z— 26-57. 



LS scheme 



Z=2Q Z=2% 



Z=32 



Z=36 



Z=42 



Z=47 



Z=54 



Z=57 



jj scheme 



3s 2 Si 



3/2 
/2 



/2 

1 S)2p 2 P 
1 S)2p 2 P 1 
3 P)2p 
3 P)2p 
3 P)2p 
3 P)2p 
1 P)2p 
1 P)2p 
3 P)2p 
1 P)2p 
3 P)2p 
3 P)2p 
3 P)2p 
3 P)2p 
3 P)2p 4 P 3/2 
1 P)2p 2 D, 
1 P)2p 2 P C 



S3/2 
4 D 5/2 
4 D 7/2 
D3/2 

2 Pl/2 

2 Dz,2 

P5/2 
Sl/2 

4 D 1/2 



/2 



l D, 



! D 5 



3/2 



12 



5/2 



3 P)2p 2 P X 
3 P)2p 2 D- 
3 P)2p 2 Si 
3 P)2p 4 P E 
1 D)2p 2 P\ 
3 P)2p 4 P 
1 D)2p 2 F, 
3 P)2p 4 P 



3/2 

/2 
3/2 
/2 
'3/2 
/2 
'5/2 
7/2 



3/2 
/2 
3/2 
'5/2 



3 P)2p 4 P X 

3 P)2 P 2 D; 

1 D)2p 2 D ; 
3 P)2p 4 D 7/2 
3 P)2p 4 D 5/2 
3 P)2p 4 D 1/2 
3 P)2p 4 S 3/2 
3 D)2p 4 P 1/2 
3 D)2p 4 P 3/2 
3 D)2p 4 F 9/2 
3 P)2p 4 D 3/2 
3 D)2p 4 P 5/2 
1 D)2p 2 F 5/2 
3 D)2p 4 F 7/2 
3 P)2p 2 S 1/2 
3 D)2p 4 F 5/2 
1 S)2p ~ 
3 D)2p 
3 D)2p 
1 D)2p 
1 D)2p 
1 D)2p 
3 P)2p 
3 D)2p 



f3l2 

4 F 3/2 
4 D r/2 
2 D 3/2 

2 Pl/2 
2f? 5/2 



'D 



5/2 



^3/2 

D)2p 4 D 1/2 
D)2p 4 D 3/2 
7/2 
5/2 
3/2 



3 D)2p 2 p 
3 D)2p 4 D, 
3 D)2p 4 F- C 
1 D)2p 2 D< 
4 S)2p 2 Pi 



/2 
/2 



-3946654 
1808945 
1909949 
2008738 
2033074 
2039466 
2039528 
2050246 
2064245 
2064483 
2081011 
2129621 
2135848 
2140088 
2148356 
2151490 
2160829 
2165759 
2234248 
2253626 
2299453 
2323992 
2326176 
2331859 
2340430 
2343264 
2360263 
2360688 
2361387 
2366738 
2367125 
2412003 
2419891 
2423029 
2432332 
2441515 
2442028 
2444121 
2446690 
2448304 
2450177 
2457171 
2458536 
2467888 
2473728 
2474319 
2475520 
2476256 
2478055 
2494618 
2503765 
2532307 
2536607 
2553254 
2553625 
2554863 
2560337 



-5408870 
2258991 
2425440 
2500357 
2527745 
2544056 
2535211 
2548295 
2571269 
2573798 
2593540 
2676539 
2685223 
2685129 
2668661 
2701237 
2733488 
2724616 
2803037 
2839597 
2882506 
2875690 
2884880 
2890891 
2900161 
2903098 
2934357 
2923345 
2925210 
2941047 
2941219 
2993290 
3004313 
2998759 
3009247 
3018138 
3044141 
3022154 
3035904 
3025385 
3066099 
3059907 
3050488 
3072651 
3057201 
3087038 
3073042 
3072416 
3117666 
3110698 
3111834 
3146832 
3128717 
3150273 
3193957 
3199650 
3171081 



-7101306 
2745096 
3005011 
3029307 
3060001 
3092521 
3068422 
3083346 
3121744 
3127096 
3150111 
3264927 
3305872 
3252230 
3231730 
3326210 
3378432 
3366873 
3437436 
3496708 
3533957 
3466402 
3488911 
3494970 
3504745 
3506655 
3558208 
3530947 
3533627 
3566519 
3566447 
3608811 
3622761 
3630191 
3646026 
3640285 
3676510 
3645035 
3658924 
3646951 
3702033 
3688115 
3688938 
3726831 
3686257 
3769447 
3769007 
3770426 
3792216 
3801760 
3774944 
3830007 
3766505 
3837776 
3905127 
3916483 
3838660 



-9720897 
3441058 
3882164 
3783964 
3820060 
3887772 
3830075 
3846505 
3919367 
3928976 
3956604 
4098156 
4239177 
4070053 
4045552 
4266154 
4358921 
4344602 
4387010 
4490270 
4517384 
4305185 
4360343 
4366156 
4375109 
4376509 
4462832 
4407111 
4410325 
4478642 
4478544 
4494602 
4511681 
4541409 
4568623 
4537586 
4576282 
4540867 
4555153 
4541206 
4615334 
4592423 
4624648 
4658307 
4592770 
4742709 
4724883 
4716884 
4823275 
4822655 
4819276 
4918354 
4684639 
4932054 
4973425 
4994914 
4841463 



-14444875 
4560251 
5436825 
4993448 
5040067 
5202696 
5053095 
5069726 
5237983 
5254252 
5289060 
5504613 
5881367 
5417213 
5387153 
5921839 
6115991 
6098245 
6072516 
6270868 
6280342 
5644107 
5793880 
5799124 
5807350 
5807902 
5937498 
5848142 
5851442 
5996208 
6004427 
5959608 
5970412 
6064401 
6048355 
6017419 
6090589 
6009260 
6029241 
6018858 
6110247 
6082264 
6178953 
6214029 
6086770 
6305020 
6270683 
6243639 
6692492 
6690006 
6600142 
6841218 
6205225 
6886124 
6875056 
6912385 
6685037 



-19131372 
5526519 
6970009 
6037431 
6094393 
6391332 
6110756 
6126870 
6430207 
6451971 
6493231 
6756597 
7491636 
6632701 
6598225 
7543576 
7878203 
7861912 
7718801 
7858384 
7824403 
6799972 
7086329 
7090560 
7097923 
7098885 
7244558 
7146773 
7149759 
7322997 
7321838 
7285858 
7282549 
7447505 
7371746 
7358627 
7459516 
7354901 
7401210 
7410665 
7495342 
7445653 
7590274 
7658925 
7440067 
7764437 
7669065 
7642346 
8552527 
8549142 
8335220 
8727106 
7605894 
8784045 
8755816 
8834753 
8542196 



-26886482 
6849732 
9000854 
7475330 
7547744 
8153954 
7569587 
7585055 
8198752 
8228075 
8278916 
8605934 
9880164 
8434674 
8393917 
9927031 
9971617 
10139433 
10074749 
10205550 
10160693 
8402456 
9174130 
9005390 
9010959 
9013658 
9239436 
9072189 
9074248 
9275443 
9262357 
9436271 
9224215 
9514794 
9328624 
9367047 
9498136 
9320567 
9431150 
9465885 
9669557 
9536164 
9666579 
9924797 
9484093 
10055141 
9823175 
9827384 
11702449 
11693655 
11178869 
11696997 
9808688 
11721257 
11862792 
11836413 
11682890 



-30654172 
7377220 
9844722 
8054553 
8133682 
8929008 
8158101 
8173589 
8976653 
9009218 
9064324 
9416474 
10785869 
9227490 
9183909 
10836191 
10881891 
11066015 
10996341 
11168436 
11394844 
9055318 
10022351 
9849121 
9853573 
9857413 
10097513 
9920622 
9921995 
10133874 
10114485 
10311330 
10078225 
10579279 
10189264 
10262012 
10404492 
10184749 
10330672 
10368561 
10796904 
10459319 
10585200 
10978300 
10396064 
11119321 
10854112 
10855818 
12906145 
12876722 
12561626 
13061250 
10842983 
13062930 
13125630 
13135682 
12867258 



3*1/2 

3si/ 2 3s 1/2 [0]2p 3/2 
3si/ 2 3s 1/2 [0]2pi /2 
3si/ 2 3pi /2 [0]2p 3/2 
3s 1 / 2 3pi /2 [l]2p 3/2 
3si /2 3p 3/2 [2]2p 3/2 
3si/ 2 3pi /2 [l]2p 3/2 
3si/ 2 3pi /2 [l]2p 3/2 
3s 1/2 3p 3 / 2 [l]2p 3/2 
3si/ 2 3p 3/2 [l]2p 3/2 
3si /2 3p 3/2 [l]2p 3/2 
3si/ 2 3p 3/2 [2]2p 3/2 
3s 1/2 3pi /2 [0]2p 1/2 
3si/ 2 3p 3 / 2 [2]2p 3/2 
3si/ 2 3p 3/2 [2]2p 3/2 
3si/ 2 3pi/ 2 [l]2pi/ 2 
3si/ 2 3p 3/2 [2]2p 1/2 
3si/ 2 3p 3/2 [l]2pi /2 
3si/ 2 3pi/ 2 [l]2pi/2 
3si/ 2 3p 3/2 [2]2pi /2 
3si/ 2 3p 3/2 [l]2pi /2 
3pi/ 2 3p 1/2 [0]2p 3/2 
3si/ 2 3d 3 / 2 [l]2p 3/2 
3pi/23p 3 / 2 [l]2p 3/2 
3f>i/23p3/2[2]2p 3/2 
3pi/23p 3 / 2 [l]2p 3/2 
3pi/23p 3 / 2 [l]2p 3/2 
3pi/23p 3 / 2 [2]2p 3/2 
3pi/23p 3 / 2 [2]2p 3/2 
3p3/2 3p 3 / 2 [2]2p 3/2 
3p3/2 3p 3/2 [2]2p 3/2 
3pi/23p 3 / 2 [2]2p 3/2 
3p3/2 3p 3 / 2 [2]2p 3/2 
3f>3/2 3p 3 / 2 [2]2p 3/2 

3p3/2 3p 3 /2[0]2p 3/2 

3si/ 2 3d5 /2 [3]2p 3/2 
3s 1/2 3d 3/2 [l]2p 3/2 
3s 1/2 3d 3/2 [l]2p 3/2 
3si/ 2 3d 3/2 [2]2p 3/2 
3si/ 2 3d 3/2 [2]2p 3/2 
3pi/ 2 3pi/ 2 [0]2p 1/2 
3s 1/2 3d 5/2 [2]2p 3/2 
3si/ 2 3d 3/2 [2]2p 3/2 
3s 1/2 3d 5/2 [2]2p 3/2 
3si/ 2 3d 5 / 2 [2]2p 3/2 
3si/ 2 3d 5 / 2 [3]2p 3/2 
3s 1/2 3d 3/2 [2]2p 3/2 
3si/ 2 3d 5/2 [3]2p 3/2 
3pi/23p 3 / 2 [2]2p 1/2 
3pi/23p 3 / 2 [l]2pi /2 
3s 1 / 2 3d 5/ / 2 [2]2p 3/2 
3pi/23p 3 / 2 [2]2p 1/2 
3s 1/2 3d 5/2 [3]2p 3/2 
3p3/23p 3 / 2 [2]2p 1/2 
3p3/23p 3 / 2 [2]2p 1/2 
3si/ 2 3d 3 / 2 [2]2p 1/2 
3pi/23p 3 / 2 [l]2pi /2 




FIG. 5: Energy splittings (AE/(Z 



— 3) 4 in cm x ) as functions of Z. 
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TABLE II: Auger energies in Fe 15+ in eV. Comparison with 
predicted data (Eat) from Bliman et al. 



LS scheme 



£rmbpt 



E, 



fit 



1 S)2p 2 P 
1 S)2p 2 P 1 



jj scheme 



3/2 
/2 



3/2 



3 P)2p 4 S : 
3 P)2p 4 D 5/2 
3 P)2p 4 D 7/2 
3 P)2p 4 D 3/2 
1 P)2p 2 P 1/2 

3/2 
Sl/2 

4 D 1/2 

4 Pl/2 
E>3/2 

2 D 5/2 

2 D 5/2 
P3/2 

2 Pl/2 
P>3/2 

2 Su 



1 P)2 P 2 D; 

3 P)2p 4 P 5 
1 P)2p 
3 P)2p 
3 P)2p 
3 P)2p 
3 P)2p 
3 P)2p 
1 P)2p 
1 P)2p 
3 P)2p 
3 P)2p 
3 P)2p 
3 P)2p 
1 D)2p 



1/2 

J 3/2 



'Pi 



3 P)2p 4 P t 
1 D)2p 2 P 
3 P)2p 4 P 



/2 
'5/2 
7/2 
3/2 



3 P)2p 4 P 1/2 

3 P)2p 2 D 3/2 

1 D)2p 2 D 5/2 
3 P)2p 4 D T/2 
3 P)2p 4 D 5/2 
P)2p 4 D 1/2 
P)2p 4 S 3/2 
D)2p 4 P 1/2 
3/2 

9/2 
3/2 
5/2 

1 J D)2p 2 P 5/2 
3 P)2p 4 P 7/2 
3 P)2p 2 S 1/2 
5/2 



£)2p 4 P 3 , 
D)2p 4 P 9 , 
P)2p 4 P 3 
! P)2p 4 P S 



£>)2p 4 P 5 



3/2 

*3/2 



'Pi 



7/2 
3/2 



/2 
5/2 



1 S)2p 
3 D)2p 
3 D)2p 
1 D)2p 
l D)2p 
1 D)2p 
P)2p 
D)2p 
D)2p 
D)2p 
D)2p 
D)2p 
D)2p 
1 D)2p 
1 S)2p 2 P 1/2 
3 P)2p 2 P 1 



2 D 5/2 
2 P 3 /2 

4 Pl/2 

4 As 



3/2 
7/2 



5/2 



3/2 



-D 



5/2 



/2 



224.3 
236.8 
249.1 
252.1 
252.9 
252.9 
254.2 
255.9 
256.0 
258.0 
264.0 
264.8 
265.3 
266.4 
266.8 
267.9 
268.5 
277.0 
279.4 
285.1 
288.1 
288.4 
289.1 
290.2 
290.5 
292.6 
292.7 
292.8 
293.4 
293.5 
299.1 
300.0 
300.4 
301.6 
302.7 
302.8 
303.0 
303.4 
303.6 
303.8 
304.7 
304.8 
306.0 
306.7 
306.8 
306.9 
307.0 
307.2 
309.3 
310.4 
314.0 
314.5 
316.6 
316.6 
316.8 
317.4 
317.5 



225.0 
237.3 
249.5 
252.4 
253.2 
253.2 
254.6 
256.0 
256.1 
258.3 
264.5 
265.1 
265.7 
266.5 
267.1 
268.1 
268.5 
277.1 
279.6 
285.1 
289.2 
289.6 
290.0 
291.3 



293.3 

294.1 
293.6 

300.4 



303.6 
303.7 
304.4 
305.3 



317.1 



3si/ 2 3s 1/2 
3s 1/2 3si/ 2 
3si/ 2 3pi/ 2 
3si /2 3pi/ 2 
3si /2 3p 3 /2 
3s 1/2 3pi/ 2 
3si /2 3pi/ 2 

3Sl/23p 3 /2 

3si /2 3p 3 /2 

3Sl/23p 3 /2 

3si /2 3p 3 /2 
3s 1/2 3pi/ 2 
3si /2 3p 3 /2 

3Sl/23p 3 /2 

3si /2 3pi/ 2 

3Sl/23p 3 /2 

3si /2 3p 3 /2 
3s 1/2 3pi/ 2 
3si /2 3p 3 /2 

3Sl/23p 3 /2 

3pi/ 2 3pi/ 2 
3si/ 2 3d 3/ /2 

3Pl/23p 3 /2 

3pi/ 2 3p 3 / 2 

3Pl/23p 3 /2 

3pi/ 2 3p 3 / 2 

3Pl/23p 3 /2 

3pi/23p 3 /2 

3P3/23P3/2 
3P3/23p 3 /2 

3pi/ 2 3p 3 / 2 

3P3/23p 3 /2 
3P3/23P3/2 
3P3/23p 3 /2 

3si/23d 5 / 2 
3s 1 / 2 3d 3/ / 2 
3si/ 2 3d 3 / 2 
3si/ 2 3d 3 / 2 
3si/ 2 3d 3 / 2 
3pi/23pi/ 2 
3si/ 2 3d5/2 
3si/ 2 3d 3; / 2 
3si/ 2 3d5/2 
3si/ 2 3d 5 / 2 
3si/ 2 3d5/ 2 
3si/ 2 3d 3/ / 2 
3si/ 2 3d5/2 
3pi/ 2 3p 3 / 2 

3Pl/23p 3 /2 

3si/ 2 3d 5 / 2 

3Pl/23p 3 /2 

3si/ 2 3d 5 / 2 

3P3/23P3/2 
3P3/23P3/2 

'3/2 

2 
2 



3si/ 2 3d : 
3pi/ 2 3p 3 / 
3p 3 /23p 3 / 



i ]2p 3 /2 

']2pi/2 
']2p 3 /2 
]2p 3 /2 
]2p 3 /2 
]2p 3 /2 
]2p 3 /2 
]2p 3 /2 
]2p3/2 
]2p3/2 
]2p 3 /2 

]2pi/ 2 

]2p 3 /2 
]2p 3 /2 

]2pi/ 2 

]2pi/2 

]2pi/ 2 
]2pi/ 2 

]2pi/2 
]2pi/2 
]2p 3 /2 
]2p 3 /2 
]2p 3 /2 

:]2p 3 /2 

]2p 3 /2 

]2p 3 /2 

]2p 3 /2 

:]2p 3 / 2 

]2p 3 /2 

:]2p 3 /2 

]2p 3 /2 

:]2p 3 /2 

]2p 3 /2 
1 ]2p 3 /2 

]2p 3 /2 

]2p 3 /2 

]2p 3 /2 
]2p 3 /2 

]2p 3 /2 

! ]2pi/ 2 
]2p 3 /2 
]2p 3 /2 

]2p 3 /2 
]2p 3 /2 

]2p 3 /2 
]2p 3 /2 

]2p 3 /2 

:]2pi/ 2 
]2pi/ 2 

]2p 3 /2 

]2pi/ 2 

]2p3/2 

]2pi/ 2 
']2pi/ 2 
]2pi/ 2 
]2pi/ 2 

»]2pi/2 



TABLE III: Wavelengths A in A for Na-like ions for odd- 
parity states given relative to the ground state. Comparison 
with experimental results from Ref. |2lj] (A a ) and Ref. |2^] 

(A*). 







Z= 


47 


Z=57 


Oonnguration 


T 

.1 


Armbpt 


A ( , 


Armbpt Af, 


3s 1/2 3s 1/2 [0]2p3 / 2 


1.5 


4.0555 


4.0547 


2.6294 2.6294 


3s 1/2 3si/2[0]2p 1/ 2 


0.5 


3.8312 


3.8309 


2.4692 


3si/23d 3 /2[2]2p 3/ 2 


1.5 


3.7423 


3.7427 


2.4249 


3s 1/2 3d 5/ 2[2]2p 1/ 2 


2.5 


3.5707 


3.571 


2.2788 


3P3/23p 3 /2[0]2p 1/2 


0.5 


3.5701 


3.571 


2.2798 


3s 1 /23d 3 /2[2]2p 1/2 


1.5 


3.5700 


3.571 


2.2766 


3Sl/23p 3 /2[l]2Sl/2 


1.5 


3.4261 


3.426 


2.1805 


3Sl/23p 3 /2[l]2Si/ 2 


0.5 


3.4121 


3.412 


2.1739 



TABLE IV: Wavelengths A in A for the 3s3p( 3 P)2p 4 D 7 / 2 - 
3s3d( 3 D)2p 4 F 9/2 transition in Na-like ions with Z=15-30. 
Comparison with experimental (A cxp t) 
data from Jupen et al. jl^] . 



and predicted (As 



z 


Armbpt 




Afit 


15 


729.74 




728.3 


16 


613.99 


614.06±0.08 


614.20 


17 


534.12 


532.70±0.10 


532.67 


18 


472.79 


471.34±0.05 


471.22 


19 


424.23 




424.05 


20 


384.80 




384.15 


21 


350.86 




351.98 


22 


325.76 


324.87±0.02 


324.89 


23 


302.61 




301.69 


24 


281.75 


281.67±0.08 


281.60 


25 


264.32 




263.98 


26 


248.73 


248.36±0.05 


248.39 


27 


234.76 




234.48 


28 


222.26 


222.00±0.10 


221.98 


29 


210.93 


210.70±0.05 


210.68 


30 


200.63 




200.39 



no experimental data between Z=47 and Z=100 for the 
sodium isoelectronic sequence. Availability of such data 
could lead to an improved understanding of the relative 
importance of different contributions to the energies of 
highly-charged ions. These calculations are presented as 
a theoretical benchmark for comparison with experiment 
and theory. The present calculations could be further im- 
proved by including third-order correlation corrections. 
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APPENDIX A: THE PARTICLE-PARTICLE-HOLE DIAGRAM 

The model space state vector for an ion with two valence electrons (v,w) and one hole (a) can be represented as 



*(QJM) = N(Q) }^(vw\K 12 }(K 12 a\K}alala a \0} 



t„t. 



(Al) 



where |0) is the state vector for the core (ls 2 2s 2 2p 6 in our case), Q describes a particle-particle-hole state with quantum 
numbers n®K®n®, K®,n ®K®, and intermediate momentum J\ 2 . We use the notation K% = {Ji,M{\ and v = {j v ,m v }. 
The sum in Eq.( [Al| ) is over magnetic quantum numbers m Vl m w , m a , and M\ 2 . The quantity (K\K 2 \K 3 ) is a 
Clebsch-Gordan coefficient: 



(K X K 2 \K 3 ) = (-1 



J 1 -J 2 + M 3 



J\ Jl J3 
Mi M 2 -M 3 



(A2) 



The above representation of the state vector is somewhat inconvenient; for example, it leads to an expression 
containing 36 terms for the second-order diagram with six free ends. It was shown in Ref. 23 that it is more efficient 
to express the state vector in a manifestly symmetric form in the case of three-electron system. Here, we can only use 

symmetry of two particles v and w. As a result, we obtain the expression for the particle-particle-hole contribution 

(2) 

E3 in second order consisting of a sum of E a , Eh, E Cl and Ed terms: 



E. 



(2) 



E„ + Eh + E r + Ed 



E a (v°w [J 12 }aJ, x y°[Ji3}bJ) = - ]T ]T P Jia (v°v, w°w)P Jl3 (xV, y°w') 



, 1 . .„ . „. , , ,, „ . .„ , . ,„, v-^ q(vwan)q(nbv'w') , . „. 

(vwK 12 (K 12 aK)(v'w'K 13 (K 13 bK) £ ^ m > , A3 

^— ' e(v) + e(w) — e(a) — e(n) 



E b (v w [J 12 }aJ,x Q y°[J 13 }bJ) = J2J2 P ^ v,w w)Pj 13 (x°v',y Q w') 



(vw\K 12 )(K 12 a\K)(v'w'\K 13 )(K 13 b\K)Y, 



g(vww'n)g(nbv'a) — g(vww'n)g(nbav') 
e(v) + e(w) — e(w') — e(n) 



(A4) 



E c (v°w [J 12 \aJ,x y [J 13 ]bJ) =^^P Jl3 (.° 1 ;,^)P Jl3 (x%' ) yV) 



vw v'w' 



{vw\K l2 ){K l2 a\K){v'w'\K l3 ){K l3 b\K)Y, 



g(vban)g(nwv'w') — g(bvan)g(nwv'w ! ) 



e(v) + e(b) — s(a) — e{n) 



, (A5) 



E d (v°w a [J 12 }a J, x°y°[ J 13 ] bJ) = ]T ]T P Jl2 (v°v, w°w)P Jl3 (xV, y°w') 

vw v'w' 

= ^^(vw\K 13 )(K i3 a\K)(v'v/\K 1 3)(K 1 3b\K)^2 



vw v'w' 



e(w) + e{b) — e(v') — e(n) 



x [— g(wbv'n)g(nvw'a) + g(bwv'n)g(nvw'a) + g(wbv'n)g(nvaw') — g(bwv'n)g(nvaw')] . (A6) 
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The quantity Pj(v°v, w°w) is a symmetry coefficient defined by 

Pj(v°v,w°w) = r] v o w o [S v o v S w o w + (-l) Jv+: ' u ' +J+1 S v o w d w o v ] , (A7) 
where rj vw is a normalization factor given by 

1 for w =/= v 
I/a/2 for W = V 

and e(i) is the one-body DF energy of t he state i. The quantities {vw\Ki%) are Clebsch-Gordan coefficients and 



(Ki2a\K) differs from the definition (A2) by the (— l)^~ m " factor and — m a in the 3j coefficient. The indices v, 
w, x, and y distinguish the valence states, indices a and b refer to the core states and index n refers to excited 
states. Two-particle matrix element g(1234) is the sum of the two-particle Coulomb and Breit matrix elements and 
is presented as a product of radial and angular parts: 

ff (1234) =^(X fe (1234) + ^(1234))(-l)i^+-^x^(-ir( ? j i M ( ^ _M . (A8) 

Radial part Xk(abcd) of the Coulomb matrix element is given by 

X k (abcd) = (a\\C k \\c)(b\\C k \\d)R k (abcd). (A9) 

The quantities Ck are normalized spherical harmonics and Rk{abcd) are Slater integrals. The expression for the Breit 
matrix elements -Bfc(1234) is given in Ref. fl28fl . Carrying out angular reduction, we obtain the final expression for the 
four contributions E a , E^, E c , and Ed to second-order diagram of particle-particle-hole interaction: 



vw v'w' 



s(v) +e(w) -e(o) - e(n) v ,JnJ \k j w j v } \ k' j w > j v < j y 



n k k 



E b {v w a [Ji2\aJ,x°y [MbJ) = £ £ Pj 12 (v°v, w°w)P Jl3 (x%\y%')y/[J^\{M 



xEEE' 

n k k' 

+Xk(nbav' 



Xy (vww'n) 



jw' jn Jl 



(v) 


+ e(w) - 


e(w') — 


e(n 


) I 3w 


jb 


jv' k 




jn 


k 1 


Jw' 


J Jl3 


} { jw' 


J 


Jl2 J 



(-1) 



i v / i2+ J i3+i a +j t ,'+i, 1 +j ll 



(-1) 



-Xk(nbv'a) 





ja 




/ jn 




k 


\ Jl2 


Jl3 




1 Jl3 


Jl2 


3 vi 



r_-\\h+jm'+jv,-J 



, (All) 



E c (v°w° [J 12 ]aJ, x°y°[J 13 ]bJ) = ^ ]T P ^ (v°v, w°w)P Jl3 (zV, y V) ^ Ji 2 ] [ Ji 3 ] 



EEE 



Xfc/ (v'w'vn) 



n k h 



jv jn J] 

-f e(v') + e(w') - e(v) - s(n) \ jw' jv' k' 
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X k (wban) 



ja jw k 1 J jb j n k 

jv J J 12 J 1 jv J Jl3 

+Xk(bwan) 



jb 


ja 


n 


/ jn 


JW 


\ 


Jl2 


Jl3 




\ Jl2 


Jl3 


Jv J 



■l\l+fc+*' 

(_X)J''+ira'+J" + -' + ' 7 12 



(A12) 
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E d (v°w [Ji2}aJ,x y°[Ji3}bJ) = E E Pj^,w°w)Pj 13 {x v' ,y W 'U[M[J^] 

vw v'w' 

- \~L+3n+jw+3a+j v <+Jl2 + Jl3 



EE' 



n kk 1 + e(b) ~ e(v') - e(n) 
Xk (nvw'a)Xk> (bwv'n) 



ja jv k 
jw J J\2 



jb jv' k' 
jw' J J13 



jw jn k 

jw' J k 



ti '«" rf *'* , «i { i ! i L } ew { £ ^:U{wS}{£ £ * } 

h jv' k' 1 , f j a i„ K If ju,' J„ K 1 f ja j„ K 1 

J Jl3 J ^ \ ju> J ^12 } \ jw J k' J \ j„ A: J 



+Xfc {nvaw')Xk' (bwv'n) 



+ x fe( nw)x fe , ( ^n) |ti ; ;, 2 }{ t y 1 } aw { t £ i } { ii 1 1> \{t 3 j : 



(A13) 
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